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Abstract 

In this paper we consider the age-related aherations of heart rate variabihty on the basis of 
the study of non-Markovian effects. The age dynamics of relaxation processes is quantitatively 
described by means of local relaxation parameters, calculated by the specific localization procedure. 
We offer a quantitative informational measure of non-Mar kovity to evaluate the change of statistical 
effects of memory. Local relaxation parameters for young and elderly people differ by 3.3 times, 
and quantitative measures of non-Markovity differ by 4.2 times. The comparison of quantitative 
parameters allows to draw conclusions about the reduction of relaxation rate with ageing and the 
higher degree of the Markovity of heart rate variability of elderly people. 
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I. INTRODUCTION 



The ageing of a human organism has been in the focus of attention in physics of hve 
systems over the past years. One of the most significant systems of vital activity of a human 
body is its cardiovascular system. Today there are a number of scientific studies on the 
problems of biological ageing of the cardiovascular human system. The latter is extremely 
sensitive to age-related as well as pathological changes in a human organism. Thus not 
only physiologists, biologists and physicians have been involved in this studies but also 
experts from other natural-science areas. Conditions of a human heart are estimated by 
means of various parameters. Thus heart rate variability (HRV) represents one of the most 
frequently used parameters, in cardiology. Nowadays there are different methods of studying 
heart rate variability dynamics. In recent years fluctuations of heartbeat dynamics have 
been studied by means of several methods derived from nonlinear dynamics and statistical 
)hysics, such as detrended fluctuation analysis (DFA) , spectral analysis |^, entropy 
^^y^l^; correlation dimension j^. In paper P| authors illustrate the problems related to 
the physiological signal analysis with representative examples of human heartbeat dynamics 
under healthy and pathological conditions which is based on two methods: power spectrum 
and detrended fluctuation analysis. In this paper different characteristics of heartbeat: 
1/f fluctuations, long-range anticorrelations (monofractal analysis), self-similar cascades, 
multifractality and nonlinearity are considered. By means of a wavelet-based multifractal 
formalism it is shown that healthy human heartbeat dynamics exhibits higher complexity 
which is characterized by a broad multifractal spectrum. In paper Q| multiresolution 
wavelet analysis has been used to study the heart rate variability in a patient with different 
pathological conditions. Noise effects of abnormal heartbeats were considered in paper . 
The correlation exceptions of heartbeat dynamics of different sleep stages often have been 



researched lately. In papers 
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13[ correlation properties of the magnitude and the sign of 



increments in the time intervals between successive heartbeats during a light sleep, a deep 
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sleep, a rapid eye movement sleep were discovered by means of the detrended fluctuation 
analysis. Multiscaled randomness Q], multifractal analysis simulation by non-linear 
oscillators fractal approach based on scaling of a frequency spectrum on power law l/cj" 
quantitative analysis are also used to analyze heart rate variability. The change 
of correlations and statistical memory effects is one of the must important questions 3] in 
heart rate variability dynamics observed with ageing 0,Q|. 

Among existing methods of researching HRV one can differentiate the methods of esti- 
mating HRV in a time area, spectral methods of estimating HRV in a frequency area, as 
well as nonlinear methods. The last group of methods has proved to be a powerful means to 
study various complex systems and has brought about significant achievements in processing 
biological and medical data. In recent years universal methods of statistical physics have 
been more often used in medicine and biology. The methods of statistical physics which 
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in the field of 



have been used to research real complex systems 
cardiology reveal essentially new opportunities for the analysis, diagnostics and forecasting 
the processes of biological ageing and diseases of a human heart. They disclose dynamic 
features of HRV, latent for classical medical methods of research. 
In this paper we offer a new method of study of the problems o 



heart activity, based on our theory of discrete non-Markov processes 



already found practical application in cardiology 



ageing of a human 
2011. This theory has 



2l|, neurophysiology [2^ 



locomotor and sensomotor activity j22], epidemiology 



the study of 



24| and seismology 



II. BASIC CONCEPTS AND DEFINITION OF THE STATISTICAL THEORY 
OF NONSTATIONARY DISCRETE NON-MARKOV PROCESSES IN COMPLEX 
SYSTEMS 

The obtained data were processed by means of the above declared technique. We use the 
results of our recent theory of discrete non-Markov random processes for the quantitative 
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description of Markovian and non-Markovian components in stochastic alteration of the 
registered data. The set of three memory functions was calculated for each sequence of the 
data. Frequency power spectra for each of these functions are obtained by using the fast 
Fourier transform. For a more detailed analysis of the properties of the system we also 
consider the frequency spectrum of the first three points of the statistical spectrum of the 
non-Markovity parameter. The spectrum of the non-Markovity parameter was introduced 



earlier in the following articles [21 
the non-Markovity parameter: 

1/2 
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In this study we use the frequency spectrum of 



dtMi{t) cos{ujt) 



cosfcjt,- 



j=0 



here i = 1,2,3... is the number of the relaxation level, is the Fourier-transform and 

a power spectrum of the ith level memory function Mi{t) (see, Eq. below). The pa- 
rameters ei{u) allow to receive the quantitative estimation of long-term memory effects in 
the experimental time series of the data as shown in Ref. From the physical point of 

view the parameter ei{uj) allows to mark out the three most important cases 3|. Markov 
and completely randomized processes correspond to values e oo, quasi-Markov processes 
(memory effects can be noticed there ) correspond to values s > 1. The limiting case with 
e ~ 1 concerns the situation with non-Markov processes, i.e., processes, where there is 
long-range memory^ 

we came to the following chain of connected non-Markov 
mr): 



In early works H 21, Q 
finite-difference kinetic equations [t 

AMn{t) 



m—1 



At 



Mnit) - rA„+i J2 M„+i(jr)M„(t - jr) 

j=0 



(1) 



Here parameters A^+i represent eigen values of the Liouville's quasioperator. The relaxation 
parameters of A„+i are determined as follows: 

.(W„_iLW„) 



An 



+1 



. (WnLW„) 



IW, 



n-1 



(2) 



The zero order memory function Mo(t) in Eq. (|T|): 



Mo(t) = a(t) = t = mT, 



Afc(O) = (5xo, Sxi, 6x2, • • • , 

-Am+fc(^) ~ {6Xm, 6Xm+l, 6Xm+2, ■ ■ ■ , 6Xm+k-l}, 

describes statistical memory in complex systems with a discrete time (A°(0) and A^_,_^(t) 
are vectors of the initial and final states of the studied system). In paper [2^ we have 
received the recurrent formula on the basis of Gram-Schmidt orthogonalization procedure, 
in which the senior dynamic variable = Wn(t) is connected with the junior one in the 
following way: 

Wo = A°(0), Wi = {zL-Ai}Wo,... (3) 

W„ = {tL - A„_i}W„„i + A„„iW„_2 + n > 1. 

The initial time correlation function (TCP) a(t) and the set of discrete memory functions 
Mn{t) in Eq. (^) are important for further consideration. The first three equations of this 
chain [t = mr is a discrete time) can be presented as follows: 

Aa(t) 



At 



- rAi ^ Mi{jT)a{t - jr) + Xia{t), 



j=Q 
m—1 



At 

j=0 



=- rA-sJ2 M,{jT)M,{t - jr) + X.M^it). (4) 

This system of finite-difference Eqs. (0} is a discrete analogue of the well-known chain 
of kinetic Zwanzig'-Mori's equations. The latter plays the fundamental role in modern 



statistical physics of non-equilibrium phenomena with a continuous time. It is necessary 
to note that the chain of Zwanzig'-Mori's equations is valid only for quantum and classical 
Hamiltonian systems with a continuous time. The chain of finite-difference kinetic Eqs. (P), 
(jU is valid for complex systems, in which there is no Hamiltonian, the time is discrete, 
and there are no exact equations of motion. However, "dynamics" and "motion" in real 
complex systems undoubtedly exist and can be registered in the experiment. The first three 
equations in the chain (j3)) form a basis for the quasihydro dynamic description of stochastic 
discrete processes in complex systems. The application of Eq. (jT)) opens up new possibilities 
in the detailed analysis of the statistical properties of correlations in complex systems. The 
existence of finite-difference Eqs. (^, (j3]) allows to evaluate unknown memory functions 
(similarly time correlation functions) directly from the experimental data. 
Let's determine the experimental relaxation time te by the equation: 



Using the experimental data now we can define the relaxation time of the studied system. 
Further we can compare the experimental time te with the theoretical one. The theoretical 
relaxation time te (where i = 1,2,3... is the number of approximation) can be determined 
on the basis of Zwanzig'-Mori's equations for various correlation approximations. For the 
first age group (young people) better accordance between the experimental and theoretical 
times is gained in the 6th correlation approximation: M^{t) = Mi{t) (see Table 1). This one 
shows the presence of long-range memory in the considered system. For the examined group 
of elderly persons better accordance of relaxation times is received in the first correlation 
approximation Mi(t) = a{t) (see Table 1). It indicates to the existence of short-range 
memory in this group. 

Using the Laplace transform on the first three Zwanzig'-Mori's equations, we shall receive: 



N 




(5) 
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sa{s) — 1 = Aia(s) — Aia(s)Mi(s), 
sMi(s) - 1 = A2Mi(s) - A2Mi(s)M2(s), 

sM2{s) - 1 = X^M^is) - A3M2(s)M3(s). 

One can solve this system by means of various approximations. For the approximation 
Mi(t) = a{t) (the first approximation in Table 1) we shall receive: 



, -{s - Ai) + ^(s-Ai)2 + 4A2 

als) = . 

^ ^ 2Ai 



Ai + ^Xj + 4A'i 



Tabic 1 

The experimental te and theoretical (i = 1, . . . ,6) relaxation times for different age groups 



Age 




Tl{Ml{t) = 


T2(M2{t) = 




T4{M2{t) = 


MM3{t) = 


TeiMsit) = 






a(t)) 


Ml it)) 


Ahit)) 


a{t)) 


a{t)) 


Ml it)) 


Young 


0.126 


0.987 


1.171 


1.172 


1.188 


1.219 


0.211 


Old 


0.147 


0.260 


7.262 


6.864 


0.789 


5.440 


0.045 



For the approximation M3[t) = Mi[t) (the sixth approximation in Table 1) we shall find: 

a(s) 



AiMi(s) + s - Ai 



2A3(s - A2) 



2A3(s - A2)(s - Ai) - Ai(A2 + {s - A2)(s - A3) - A3) + Aiv/(A2 + (s - A2)(s - A3) - A3)2 + 4A3(.s - A2)(s - A3) 
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. 2A2A3 

tq = limafsj 



^-0 ^' Ai(A2 + A2A3 - A3) - 2A1A2A3 - Ai V(A2 + A2A3 - A3)2 + 4A2A3A3' 
In Table 1 we show the experimental times of relaxation te (see Eq. (5)) and theoretical 
times of relaxation Tj for different age groups. 



III. EXPERIMENTAL DATA 



We used the time series of R-R intervals in young and elderly subjects as the experimental 
data Two groups of healthy human subjects: 10 young (mean age 27 yr, range 21-34 
yr) and 10 elderly (mean age 74 yr, range 68-81 yr), participated in this study. Each group 
consisted of five women and five men. All subjects provided written informed consent and 
underwent a screening history, physical examination, routine blood count and biochemical 
tests, electrocardiogram, and exercise tolerance test. Only healthy, nonsmoking subjects 
with normal exercise tolerance tests, without any medical problems, and being on no medi- 
cation were admitted to the study. 

All subjects remained in an inactive state in sinus rhythm while watching the movie 
"Fantasia" (Disney) to maintain wakefulness. Each heartbeat was annotated by means of 
an automated arrhythmia detection algorithm, and each beat annotation was verified by 
visual inspection. The R-R interval (interbeat interval) of time series for each subject was 
then computed ^. 



IV. DISCUSSION OF THE RESULTS 



The basic outcomes are submitted in this section. Further the appropriate analysis of the 
experimental data will be carried out both for young and old people. Two new qualitative 
procedures of the appropriate analysis have been used. The procedure of the window- 
time behavior shows great oscillations of an R-R interval for the power spectra of memory 
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functions with respiratory arrhythmia. The calculation of the local relaxation parameters is 
carried out by means of the procedure of time localization. The decrease of the relaxation 
rate in ageing people is indicated through the time dependence of the localized relaxation 
parameters. The quantitative assessment of non-Markovity effects of heart rate variability 
is carried out by means of a special measure. 

A. Study of age-related alterations of heart rate variability 

Further we submitted figures for one young and one old person. The figures reflect a 
general pattern of the group. Below we show the analysis which describes the experimental 
data for the first and second group. 

In Fig. 1 we have presented the time series of the first four dynamic orthogonal variables 
Wi, where i = 0, ...,3 for young (Figs, la-d) and elderly (Figs, le-h) people. The time 
series of the elderly people differ in greater amplitude and frequency of fluctuations of an 
R-R interval. The scattering interval of the oscillations of an R-R interval for young people 
constitutes (0.83r ^ 1.18r, where r = 1.0244s, r-time of discretization). The variability of 
an R-R interval in elderly people changes within the limit of (0.92r^l.5r, where r = 1.062s) 
for the initial signal. The frequency of fluctuations of an R-R interval in elderly people is 
greater than in young people. The normal cardiac activity of elderly people is accompanied 
by more randomized fluctuations of an R-R interval. 

In Fig. 2 we have presented the phase portraits for the variability of an R-R interval 
in young and elderly people on six plane projections of the first four dynamic orthogonal 
variables. The phase clouds of the young and elderly people are symmetric concerning the 
origin of coordinates frame. The phase clouds have a centered nucleus. The nucleus encloses 
some points, speckled on the perimeter. The phase points for elderly people have a 2-3 times 
greater interval of a scattering, due to the presence of more appreciable oscillations of an 
R-R interval. 
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Figure 1: The time series of the first four dynamic orthogonal variables Wi, where i = 0, ..,3 for the 
young (a-d) and the elderly persons (e-h). The time series of the elderly person is accompanied by 
more significant fluctuations of an R-R interval. The time series of elderly person is characterized 
by greater frequency of occurrence of significant fluctuations. It shows higher degree of Markovity 
of the fluctuations of an R-R interval for the elderly person. 

In Fig. 3 we have presented the power spectra of the initial TCF /io(cu) (Fig. 3a) and 
the first three memory functions of younger orders fii{uj) (where i = 1,2,3) (Figs. 3b-d) 
for young and elderly people. The figures are submitted on a log-log scale. The power 
spectrum of the initial TCF for the elderly person differs in smaller dynamic fractures in 
the area of low frequencies. It is possible to find dynamic splashes (dynamic peaks) on all 
figures in the area of frequencies (0.2/.M. < u < 0.5 f.u., where If.u. = 1/t), in particular. 
These dynamic peaks appear due to respiratory arrhythmia. The given dynamic splashes 



remind of the well known shape of the Suyumbike Tower 



2l| . The increase of the power 



spectrum on these frequencies reflects age-physiological changes. Ageing people develop a 
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Figure 2: The phase portraits of variabihty of an R-R interval on six plane projections of various 
combinations Wi,Wj for young and elderly people. The phase clouds of the young man are more 
compressed towards the center of coordinates system. Around the centralized nucleus of the phase 
clouds for the elderly person, separate points are scattered. The interval of their disorder 2-3 times 
exceeds the areas of disorder of the phase points for the young man. 

shift of these dynamic peaks in the range of high frequencies. The dynamic peaks, which 
relate to the respiratory arrhythmia of the young man, are discovered within the frequency 
interval of 0.25/.n. < uj < 0.45/.n. These dynamic splashes of the elderly person are in the 
range of OAf.u. < uj < 0.55/.n.. The specific arrhythmia of cardiac activity at respiration 
accounts for this conclusion. This frequency (at respiratory arrhythmia) in elderly persons 
is higher than in young ones. 

The procedure construction of the window-time behavior of the power spectra of memory 
functions leads to the similar conclusion. The similar procedure allows to consider in detail 
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Figure 3: The power spectra of the initial TCF (a) and the first three memory functions (b-d) for 
the young and elderly persons. The power spectra of the initial TCF of the elderly person differs in 
smaller dynamic breaks in the field of low frequencies. In all the power spectra the dynamic splashes 
(peaks) are found in the field of high frequencies, which is connected with respiratory arrhythmia. 
These dynamic peaks of the elderly person are shifted in the area of higher frequencies. It testifies 
to the increase in frequency of cardiac reductions at breath that comes with ageing. 



any dynamic regularities originating in the power spectra of the memory functions. The idea 
of this procedure consists in the following [2^. Originally it is necessary to determine the 
optimal length of the sample. When the length of the sample is small, the "accumulated" 
information will be insufficient for carrying out a qualitative correlation analysis due to 
gross errors and the influence of noise effects. When the length of the sampling is large the 
necessary "sensitivity" weakens. The analysis of samples of different lengths shows, that 
the optimal length for this procedure constitutes 2^ = 256 points. From the initial array of 
the experimental data we take 256 initial points. We receive the first window of 256 points. 
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Then we build a memory function power spectrum for this sample and take the next window 
of 256 points (from 257 to 512). Then we build the power spectrum of the memory function. 
This procedure is carried out repeatedly up to the end of the sampling of the experimental 
data. In Fig. 4 the time- window behavior of the first memory function fii{uj) in young and 
elderly people is submitted. The most appreciable dynamic splashes (peaks) in the power 
spectra are connected with respiratory arrhythmia. All dynamic peaks are found in the 
particular range of frequencies. Generally in young people the range of these dynamic peaks 
meets 0.25/.U. < uj < 0.45/.U.. For elderly people this range is shifted to the right and 
meets OAf.u. < u < 0.55/.M. This implies amplification of cardiac activity at respiration 
with ageing. 

B. Quantitative measure of the efTects of non-Markovity in heart rate variability 

Any complex system has a great number of degrees of freedom. The high dimension of 
complex systems, the presence of strong nonlinear interactions and the feedback determine 
their behavior. This behavior can be characterized by Markov random processes. Strong 
external influence at accidents, crises and human diseases entails partial synchronization of 
natural chaotic behavior of complex systems. This synchronization results in the compelled 
organization of the structure of a real system and the occurrence of regular communica- 
tions. The behavior of the system becomes more ordered. Such behavior is defined by the 
amplification of non-Markov statistical effects. 

The basic idea of our method consists in defining the quantitative proportion between 
Markov and non-Markov effects of the studied stochastic process. As a quantitative measure 
of non-Markovity we suggest using the first point of the non-Markovity parameter ei{uj), 
where u = Of.u. The physical sense of this parameter consists in distinguishing Markov (pro- 
cesses with instant or short memory) and non-Markov (processes with long-range memory) 
stochastic processes. The increase of this parameter (ei(0) >> 1) means greater Markovity 
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Figure 4: The window-time behavior of the first memory function power spectra for young 

and elderly people. On certain frequencies the dynamic peaks are distinctly visible in both power 
spectra. These peaks are connected with respiratory arrhythmia. These peaks are shifted towards 
higher frequencies for elderly people. It confirms the conclusion made as a result of the analysis of 
the previous figure. 

of conditions of the system, the reduction of this parameter describes the processes with 
statistical effects of non-Markovity. The increase or decrease of the non-Markovity param- 
eter makes it possible to judge the degree of non-Markovity in the complex system. Thus 
non-Markov effects are characterized by the presence of long-range correlations and the 
amplification of statistical displays of long-range memory. Markovian processes are charac- 
terized by the amplification of short-range correlations and the reduction of the effects of 

1 

statistical memory. The quantitative measure of non- Markovity £i(0) = {ji7§y}' reveals 
the nature of the behavior of a system. 

In Fig. 5 the frequency dependence of the first three points of the non-Markovity pa- 
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Figure 5: The frequency dependence of first three points of the non-Markovity parameter £i{u) for 
the young and elderly persons. The value of the first point of the non-Markovity parameter £i{0) 
(a) on zero frequency is an original quantitative measure of non-Markovity of the process under 
study. The value of this parameter for the young man is equal to 6.24. For the elderly person the 
value of this parameter is equal to 23.16. The comparison of these values indicates, that the heart 
activity of the elderly person is characterized by the greater Markovity. 

rameter ei{u!), where z = 1,2,3 in young and elderly people, is submitted. The value of a 
quantitative measure of the degree of non-Markovity of the young man is £:i(0) = 6.24. The 
value of this parameter of the elderly person is £i(0) = 23.16. The ratio of the quantitative 
measures of non-Markovity is 3.7 times. It testifies to the increase of Markovian effects in 
fluctuations of an R-R interval with ageing. Markovity of heart rate variability is connected 
with significant fluctuations in the initial signal of the elderly person. 

Further we have presented statistical results of processing for the first and second age 
groups. In Fig. 6 we have presented the first points of the non-Markovian parameter, 
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averaged for the groups for ten young and ten elderly people. The frequency dependence of 
such parameters is defined as follows: 

N 

10 

^ly^jav.val — — ; ••• 

The value of the first point of the non-Markovity parameter ei{0) av.vai of young people is 
6.41. The value of this parameter for elderly people is equal to 27.41. The ratio of these 
values is 4.2 times. It testifies to the amplification of Markov effects in the heart activity 
in elderly people. Therefore the number of Markov components affecting the heart activity 
increases with ageing. 



C. Age-related alterations of relaxation modes 

The local relaxation parameters allow to estimate the relaxation rate in the systems. 
The procedure of localization enables to reveal internal peculiarities of the dynamics of 
cardiac activity, latent for usual correlation analysis. The idea of the method consist in the 
following. From the initial time series we take a sampling N points in length for which we 
calculate kinetic and relaxation parameters. Then we performed "step-by-step shift to the 
right" operation one interval to the right and calculate kinetic and relaxation parameters. 
This procedure is carried out up to the end of the time series. Thus received local relaxation 
parameters have high sensitivity effects of alternation and non-stationarity. If there is any 
irregularity in the initial time series it will be instantly revealed in the time behavior of local 
parameters. 

When using this method (as well as in the first procedure) it is necessary to define the 
length of the sample which allows to receive the most trustworthy information. As a result 
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Figure 6: The frequency dependence of the first point of non-Markovity parameter, averaged on 
the group of young (a) and elderly (b) people. These figures allow to define the generalized degree 
of Markovity and non-Markovity for the first and second age groups. The ratio of quantitative 
measures of the degree of non-Markovity 5 = for the second and first age groups consti- 

tutes 4.2 times. This implies that the variability of an R-R interval of elderly people becomes more 
Markovian (for the whole group). The heart activity of young people becomes more ordered and 
is characterized by high regularity. 

of the research of different lengths of the local sample we have calculated the optimal length 
which contains 2^ = 128 points. 

In Fig. 7 we have presented the time dependence of local kinetic and relaxation param- 
eters Aj, where i = 1,2,3 and Ai,A2 averaged for ten young and ten elderly people. The 
physical sense of the first relaxation parameter Ai consists in defining the relaxation rate 
for the studied process. On average the amplitude of local parameter Ai for young people 
changes within the interval of 0.1719r~^ < Ai < 0.8522r~^. For elderly people this variable 
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Figure 7: The time dependencies of local kinetic and relaxation parameters Ai, A2, A3 and Ai, A2, 
averaged on the whole group for young and elderly people. The physical sense of the first relaxation 
parameter Ai consists in defining the relaxation rate. The ratio of root-mean-square amplitude of 
this parameter to the first and second age groups is equal to 3.3 times. It indicates to higher 
relaxation rate of cardiac activity for the young people. 



changes within the interval of 0.0187r ^ < Ai < 0.3580r ^. The ratio of root-mean-square 



JV-l N 2 



amplitude (A) 



E 



AT 



for young and elderly people is equal to 3.3 times. In Table 



N-l 

2 we present the root-mean-square amplitude, dispersion cr^ = ^ 5x| and a root-mean- 

( N-l 1 2 

square deviation cr = < jj ^ {xj — (x))^ > for local relaxation parameters. The comparison 

I ^=0 J 

of these characteristics indicates reduction of the relaxation rate with ageing (the increase of 
the relaxation time). For example, the difference of the amplitudes of relaxation parameters 
Ai constitutes 1.7 times. 
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Table 2 

Some kinetic and relaxation parameters (absolute values) for young and elderly people, calculated from 
our theory 



Age 


young 


old 


young 


old 


young 


old 


young 


old 


young 


old 


Parameter 




1 




1 


A3(r-^) 






) 




) 


(A) 


0.4551 


0.1384 


0.0490 


1.1021 


1.009 


1.004 


0.083 


0.048 


0.083 


0.107 




0.016 


0.003 


0.004 


0.010 


3*10^^ 


5*10^ 


0.004 


0.001 


0.004 


0.011 


a 


0.128 


0.054 


0.060 


0.101 


0.018 


0.024 


0.063 


0.036 


0.062 


0.107 



Thus, using different approaches in the research of the time series of an R-R interval we 
have arrived at the general conclusion. The work of a heart becomes more Markov (the 
effects of statistical memory disappear) with ageing and the speed of relaxation of cardiac 
activity is reduced. 



V. CONCLUSION 

The achieved results allow to come to the following conclusions. The procedure of lo- 
calization makes it possible to calculate quantitative characteristics describing the speed of 
relaxation of cardiac activity. We have revealed essential distinctions in relaxation processes 
for different age groups on the basis of the comparative analysis of quantitative characteris- 
tics of variability of an R-R interval. The processes occurring in the heart of a young person, 
have a greater relaxation rate. Hence, at the appearance of any destructions in cardiac ac- 
tivity, faster restoration of its usual normal mode can be observed. The relaxation rate of 
cardiac activity decreases in elderly people. Heart activity comes back to its normal rhythm 
slower in this case. 

The use of the first point of the non-Markovity parameter ei (u) allows to estimate quan- 
titatively Markovian and non-Markovian effects of heart rate variability. The work of a heart 
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of an older person is characterized by greater Markovity due to the influence of a greater 
number of components, that reduce the effects of long-range memory (deterioration of an 
organism, ageing and physiological changes of a heart and other life-support systems etc.). 
The comparison of values of the parameter ei{uj) for young and elderly people (their ratio 
is equal to 4.2 in case oi uj = Of.u.), indicates a high degree of Markovity of heart rate 
variability of elderly people. The heart rate variability of young people is characterized by 
greater non-Markovity, that indicates the smaller number of its Markov components. Thus, 
the number of Markov components of cardiac activity increases with age. 

These conclusions are interconnected and supplement each other. The increase of the 
number of Markov components, affecting cardiac activity, results in the increase of the 
relaxation time of a system. The system needs greater time for restoration to the normal 
operating mode. It is caused by the decrease of long-range correlations and reduction of the 
effects of statistical memory with ageing. On the contrary, the increase of the relaxation 
rate of a system testifies to the increase of the number of regular components. For example, 
the effects of statistical memory and long-range correlations are amplified the dynamics of 
HRV. A higher relaxation rate is characteristic of normal heart activity. 

The procedure of the window-time behavior allows to find out additional age features of 
cardiac activity of a person. The frequency of cardiac reductions at breath increases with 
age. It shows the age-related displacement of dynamic bursts (connected with respiratory 
arrhjd;hmia) in the area of higher frequencies. 
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